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Preface 

This  report  is  concerned  with  heat  transfer  rates  from  a  vibrating 
cylinder  in  free  convection  in  air.  It  has  not  only  been  an  interesting 
experience  in  heat  transfer,  but  also  a  long,  laboring  experience  in 
vibration  techniques  for  obtaining  higher  and  higher  vibration  inten¬ 
sities  from  a  cylinder  when  excited  at  an  end.  The  higher  vibration 
intensities  were  possible  but  they  were  always  in  three  degrees  of 
motion,  which  for  the  purposes  of  this  investigation  was  not  applicable. 

I  wish  to  express  my  appreciation  to  those  members  of  the 
Mechanical  Engineering  Department  who  graciously  gave  their  time  and 
energy.  I  am  especially  grateful  to  Dr.  Andrew  J.  Shine  for  suggesting 
this  thesis  topic.  As  my  thesis  advisor,  Dr.  Shine  gave  timely  sugges¬ 
tions  that  were  invaluable.  Mr.  Flahive  and  Mr.  Brown  assisted  me  with 
many  practical  problems  associated  with  the  test  equipment. 

I  wish  to  thank  the  Electrical  Engineering  Department  for  the  loan 
of  equipment  and  calibration  of  the  electrical  instruments  used  in  the 
experiment . 
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Abstract 

Methods  to  increase  the  heat  transfer  rate  from  a  body  to  its 
surroundings  include  the  effects  of  vibration.  This  study  is  an  inves¬ 
tigation  of  the  effect  of  sinusoidal  vibration  in  a  vertical  plane  on 
the  heat  transfer  rate  from  a  cylinder  in  free  convection  in  air. 

Three  different  cylinders  with  diameters  0.072  in,  0.120  in,  and 
0.25  in  were  vibrated  at  amplitudes  from  0  to  1.5  in,  at  frequencies 
from  16  to  80  cps,  and  at  surface  temperatures  of  100  F  and  200  F  above 
room  temperature. 

The  heat  transfer  rate  changes  from  that  obtained  in  a  free  con¬ 
vection  process  to  that  obtained  in  a  forced  convection  process  at  the 
conventional  forced  convection  curve  of  McAdams  and  generally  parallels 
this  curve  to  Reynolds  numbers  as  high  as  1100  which  was  the  highest 
obtained.  For  a  given  Reynolds  number,  the  heat  transfer  rate  is 
independent  of  the  cylinder  temperature  and  test  position  along  the 
length  of  the  cylinder.  The  results  of  this  investigation  agree 
closely  with  results  previously  presented  for  cylinders  vibrated  in 
free  convection  in  air,  regardless  of  the  mode  or  method  of  vibration. 

The  scope  of  the  investigation  was  limited  by  vibrating  character¬ 
istics  of  the  cylinders,  failure  of  the  cylinders  at  higher  temperatures, 
and  the  limited  mass  of  the  mounting  assembly  for  the  cylinders. 
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EFFECT  OF  VIBRATION  ON  HEAT  TRANSFER  FROM 
CYLINDERS  VIBRATED  SINUSOIDALLY  WITHIN 
A  VERTICAL  PLANE  IN  FREE  CONVECTION 


I.  Introduction 


Background 

Free  convection  is  conduction  of  heat  from  a  body  to  its  boundary 
layer  with  subsequent  fluid  movement  caused  by  changes  in  fluid  density 
in  the  boundary  layer.  The  boundary  layer  is  an  effective  barrier  to 
the  transfer  of  heat  from  a  body  vo  its  surroundings;  and  to  increase 
the  rate  of  heat  transfer,  the  boundary  layer  must  be  altered.  Common 
techniques  employed  to  alter  the  boundary  layer  are  to  reduce  the  thick¬ 
ness  of  the  boundary  layer,  increase  transverse  fluid  motion  in  the 
boundary  layer,  or  both.  Vibration  is  one  method  of  altering  the 
boundary  layer  by  mechanically  inducing  turbulence  in  the  boundary  layer 
and  by  making  the  boundary  layer  thinner.  This  is  illustrated  in  the 
Schlieren  photographs  in  Figure  13. 

Within  the  past  twenty  years,  a  number  of  investigations  have  been 
performed  to  study  the  effect  of  various  methods  of  vibration  on  the 
heat  transfer  rate  from  flat  plates  and  cylinders  in  free  convection  in 
air.  The  following  conclusions  were  made  from  these  investigations: 

1.  The  heat  transfer  rate  increases  with  vibration  beyond  a 
critical  vibration  intensity. 
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2.  The  effect  of  vibration  is  independent  of  temperature  in  the 
range  of  temperature  tested. 

3.  The  heat  transfer  rate  is  independent  of  the  mode  of  vibration, 
and  is  only  dependent  on  the  intensity  of  vibration  (Ref  6:3). 

Knowledge  of  heat  transfer  rates  as  a  function  of  vibration  intensity 
has  possible  applications  in  the  design  of  equipments  which,  contain 
components  operating  at  elevated  temperatures  and  vibrating  at  known 
intensities. 

In  1962,  Shine  and  Jarvis  investigated  the  effects  of  vibration  on 
the  heat  transfer  rate  from  cylinders  using  a  different  method  of 
vibrating  the  cylinders  than  had  been  previously  used.  They  mounted  the 
test  cylinders  horizontally  at  the  ends  and  excited  the  cylinders  at  a 
point  near  one  end.  An  investigation  of  the  effect  of  vertical  vibration 
on  the  heat  transfer  rate  from  cylinders  vibrated  in  a  sinusoidal  wave 
form  in  air  was  made.  The  cylinder  diameters  were  0.032  in  and  0.072  in. 
The  frequency  range  was  from  15  to  75  cycles  per  second  and  the  amplitude 
range  was  from  0.002  to  0.99  in.  They  reported  that  the  variation  of 
the  heat  transfer  coefficient  with  vibration  intensity  generally 
paralleled  the  forced  convection  curve  recommended  by  McAdams  in  the 
range  of  vibration  intensities  tested  (Ref  8:2). 

This  investigation  was  an  extension  of  the  work  of  Shine  and  Jarvis 
to  study  the  effects  at  higher  Reynolds  numbers,  at  two  cylinder  surface 
temperatures,  and  at  different  test  positions  along  the  cylinder. 

Reynolds  number  as  used  throughout  this  report  to  express  vibration 
intensity  is  computed  from  2afdp/p,.  The  vibration  intensity  for  a 
given  cylinder  diameter  is  defined  as  2af  which  is  in  fact  the  average 
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velocity  of  the  cylinder  action.  Therefore,  the  vibration  intensity 
(2af)  and  Reynolds  number  for  a  vibrating  cylinder  are  related  because 
the  vibration  intensity  is  used  in  computing  the  Reynolds  number. 

Purpose 

The  purpose  of  this  investigation  was  to  determine  the  effect  of 
vibration  on  the  heat  transfer  rate  from  a  heated,  horizontally-mounted 
cylinder  vibrated  in  a  sinusoidal  wave  form  in  a  vertical  plane  in  free 
convection.  The  specific  objectives  of  this  investigation  were  to 

1.  Investigate  the  effect  at  Reynolds  numbers  higher  than 
previously  obtained  and  correlate  the  results  with  the 
conventional  forced  convection  curve  of  McAdams  (Ref  5:275). 

2.  Determine  the  effect  of  the  location  of  the  test  position 
along  the  cylinder  during  vibration  in  a  sinusoidal  wave  form. 

3.  Determine  the  effect  of  cylinder  temperature. 

4.  Determine  the  effect  of  cylinder  diameter. 

Scope 

This  investigation  was  concerned  with  heat  transfer  from  cylinder 
diameters  of  0.072  in,  0.120  in,  and  0.25  in. 

The  cylinders  were  vibrated  over  a  frequency  range  from  16  to  80 
cycles  per  second  and  an  amplitude  range  from  0  to  1.5  in.  The  0.072  in 
cylinder  was  tested  at  a  surface  temperature  100  F  above  room  temperature 
and  the  0.120  in  and  0.25  in  cylinders  were  tested  at  100  F  and  200  F 
above  room  temperature.  The  heat  transfer  rates  were  determined  at 
different  positions  along  the  cylinder  while  it  was  vibrated  in  a 
sinusoidal  wave  form.  Schlieren  photographs  of  the  boundary  layer  were 
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?-**£»  to  correlate  toe  behavior  of  the  boundary  layer  with  heat  transfer 
•lata  in  the  vicinity  of  the  critical  Reynolds  number* 
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The  cylinder  was  vibrated  at  a  frequency  and  amplitude  while  the 
surface  temperature  of  the  cylinder*  which  was  electrically  heated*  was 
maintained  at  a  predetermined  constant  temperature  above  roars  temperature. 
Tie  electrical  power  input  to  the  cylinder  provided  a  measure  of  the  heat 
flux  from  the  surface  of  the  cylinder,  Kussali  numbers  were  computed 
from  the  power  input  data  and  Reynolds  numbers  were  computed  from  the 
vibration  Intensity  data.  These  dimensionless  parameters  were  plotted 
to  present  the  results  in  graphical  form  for  correlation  with  previous 
investigations . 
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II,  Experimental  Equipment 

lest  Cylinders 

The  test  cylinders  were  stainless  ste^l  tubing  of  three  different 
diameters.  The  lengths  of  the  cylinders  were  also  different  because 
the  vibration  responses  differed,  and  the  lengths  were  chosen  to  obtain 
maximum  amplitude  response.  They  are  listed  below: 


Body 

Surface 

Diameter 

Length 

Wall  Thickness 

Stainless 

Steel 

Folished 

0.072  in 

O.d 

25.5  in 

0.020  in 

Stainless 

Steel 

Polished 

0,120  in 

O.d 

37.47  in 

0.020  in 

Stainless 

Steel 

Gray-oxidized 

0.250  in 

O.d 

39.17  in 

0.020  in 

A  movable,  iron-constantan  thermocouple  was  threaded  through  the 
opening  of  the  cylinder  with  one  lead  extending  out  of  each  end. 

Copper,  electrical  connectors  were  fitted  tightly  to  the  extreme  ends 
of  the  cylinders  to  provide  connections  for  electrical  power  input  and 
voltage  pick  offs  (see  Figure  3).  The  cylinders  were  mounted  at  points 
immediately  to  the  inside  of  the  connectors,  and  the  vibration  input 
shaft  from  the  vibrator  to  the  cylinder  was  fastened  to  the  cylinders 
at  a  point  6  in  from  the  left-end  mounting  (see  Figure  3). 

Mounting  Assembly 

The  mounting  assembly  for  the  test  cylinders  consisted  of  two 
vertical  posts  fastened  to  a  horizontal  beam  which  was  supported  at 
each  end  (see  Figure  3).  The  space  above  the  mounting  assembly 
including  the  test  cylinder  was  enclosed  with  a  rectangular  hood  which 
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measured  48  x  12  x  30  in.  The  bottom  and  top  of  the  hood  were  vented 
to  allow  free  convection  around  the  cylinder  while  protecting  the 
cylinder  from  random  convective  currents  in  the  room  (see  Figure  4). 

An  adjustable  weight  was  attached  to  the  right  end  of  the  test  cylinder 
through  a  pulley  arrangement  to  compensate  for  thermal  expansion  and 
thus  maintain  a  constant  tension  on  the  cylinder. 

Vibrator 

The  source  of  energy  for  vibrating  the  cylinder  was  an  electro¬ 
magnetic  exciter  having  a  variable  output  of  up  to  25  pound  force.  The 
exciter  was  powered  by  a  power  supply  which  contained  an  audio  frequency 
oscillator  and  a  variable  power  output  unit.  Frequency  of  vibration 
was  controlled  by  the  oscillator,  and  the  amplitude  of  vibration  was 
controlled  by  the  variable  power  output  rheostat.  The  exciter  was 
positioned  directly  below  the  test  cylinder,  and  was  connected  to  the 
cylinder  by  a  3/8  in  diameter,  12  in  long  aluminum  drive  rod  (see 
Figure  3).  The  rod  constrained  the  cylinder  to  vibrate  within  a  vertical 
plane  in  a  sinusoidal  wave  form. 

Heat  Power  Source 

Direct  current  electrical  power  for  heating  the  test  cylinder 
was  taken  from  a  28  v,  540  amp  rectifier  which  serves  as  the  central 
source  of  D.C.  power  for  the  M.  E.  Laboratories  at  AFIT.  The  output 
voltage  was  adjustable  from  0  to  50  volts,  and  to  provide  finer  adjust¬ 
ments  in  power  through  the  test  cylinder,  four  5.8  ohm,  10  amp,  variable 
resistors  were  connected  in  parallel  and  placed  as  a  unit  in  series  with 
the  power  circuit  to  the  test  cylinder. 
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Measurement  Equipment 

Electrical  power  to  heat  the  cylinders  was  measured  on  a  Weston 
Electric  model  931  ammeter  and  voltmeter. 

Frequency  of  vibration  of  the  cylinder  was  measured  with  a  strobotac 
model  number  631-B  manufactured  by  General  Radio  Company. 

Amplitude  of  vibration  was  measured  from  photographs  taken  of  the 
vibrating  cylinder  motion  using  a  light  source,  mirrors,  and  a  camera. 

The  arrangement  of  these  equipments  is  shown  in  the  schematic  drawing 
of  Figure  1. 

The  source  of  light  was  a  mercury  lamp  unit  manufactured  by  the 
George  W;  Gates  Company.  The  lamp  was  hooded  .and  its  light  passed 
through  a  condensing  lens  onto  a  1/32  in  pin  hole  in  the  forward  end 
of  the  hood.  This  provided  a  point  source  of  light  which  was  passed 
onto  a  series  of  mirrors,  across  the  cylinder,  and  onto  a  photographic 
film. 

Two  5.5  in  diameter,  35.5  in  focal  length,  parabolic  mirrors  and 
a  2  in  diameter  flat  mirror  mounted  on  the  rotating  shaft  of  a  small 
D.C.  electric  motor  were  used  to  direct  light  across  the  cylinder  and 
onto  the  photographic  film.  The  first  parabolic  mirror  was  positioned 
so  that  it  directed  the  light  across  the  test  cylinder  in  a  parallel 
light  beam.  The  second  parabolic  mirror  was  positioned  on  the  opposite 
side  of  the  cylinder  from  the  first  mirror.  It  received  the  light 
from  the  first  mirror  and  passed  it  onto  the  2  in  rotating  flat  mirror 
which  was  located  at  the  focal  point  of  the  second  parabolic  mirror. 

The  rotating  mirror  reflected  the  light  in  a  sweeping  motion  across  the 
aperture  to  the  camera  (see  Figure  5). 
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A  Polaroid-back  camera  containing  type  47  Polaroid  film  was  used 
to  photograph  the  motion  of  the  cylinder.  The  camera  was  mounted  on 
the  aft  end  of  a  special  plywood  housing  with  an  extendable  bellows 
section  which  contained  the  aperture  to  the  camera,  and  which  prevented 
stray  light  from  entering  the  camera. 

Temperatures  Were  measured  with  i.-on-ccnstantan  thermocouples  and 
a  Bristol's  Dynamaster  Recorder  potentiometer.  Iron  and  constantan 
wire  with  varnished  fiberglass  sheathing  was  used  in  constructing  the 
thermocouples.  The  junction  was  formed  by  spot  welding  the  ends  to¬ 
gether.  The  thermocouple  used  inside  the  test  cylinder  was  further 
sheathed  in  teflon  spaghetti  except  for  the  junction  which  protruded 
through  a  hole  in  the  spaghetti.  The  thermocouples  and  recorder  as  a 
combination  were  calibrated  against  an  accurate  mercury-in-glass 
thermometer  using  an  ice  bath,  ambient  air,  and  boiling  water  as 
references . 

Boundary  Layer 

The  boundary  layer  was  photographed  using  Schlieren  photography, 
and  the  same  equipment  used  to  photograph  the  cylinder  motion  was  used 
with  the  following  exceptions.  The  light  source  was  a  locally  constructed 
condenser  discharge  spark  unit  and  an  ionizable-gas-filled  lamp.  The 
lamp  was  housed  in  the  same  hood  used  with  the  mercury  lamp  and  was 
adjusted  within  the  hood  in  the  same  way  as  the  mercury  light.  The 
hole  in  front  of  the  light  hood  served  as  the  first  knife  edge  and  an 
adjustable  knife  edge  was  used  at  the  focal  point  of  the  second 
parabolic  mirror. 
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III.  Measurements 

The  raw  data  required  in  the  computation  cf  the  dimensionless 
parameters,  which  are  used  to  present  graphically  the  results  of  this 
investigation,  were  the  variable  quantities  of  ambient  air  temperature, 
cylinder  surface  temperature,  frequency  of  vibration,  amplitude  of 
vibration,  electrical  current  and  voltage,  and  barometric  pressure. 

The  measurement  of  these  variables  is  discussed  in  the  following 
paragraphs.  An  analysis  of  errors  in  variables  is  presented  in 
Appendix  D. 

Temperature 

The  two  temperatures — ambient  air  and  cylinder  surface — were 
measured  to  the  nearest  1/2  F  with  an  iron-constantan  thermocouple 
connected  to  an  automatic  recording  potentiometer.  The  ambient  tempera¬ 
ture  was  measured  at  a  point  level  with  the  test  cylinder  and  6  in 
away  in  a  direction  perpendicular  to  the  plane  of  vibration.  At  this 
distance  the  thermocouple  was  outside  the  effect  of  the  temperature  of 
the  cylinder.  Ambient  temperature  measurements  were  recorded  once  each 
minute  during  a  test  run. 

The  surface  temperature  at  any  test  point  along  the  cylinder  was 
measured  from  the  inside  surface  of  the  cylinder  with  a  movable  iron- 
constantan  thermocouple.  Since  the  entire  thermocouple  was  within  the 
cylinder,  it  did  not  affect  the  boundary  layer.  The  thermocouple  was 
connected  to  15  consecutive  channels  of  the  recording  potentiometer 
which  recorded  each  4  seconds  giving  an  essentially  continuous  indication 
of  the  temperature  to  assist  in  adjusting  the  power  input.  During  test 
runs,  the  ends  of  the  cylinder  were  closed  with  a  pliable  material  to 
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prevent  heat  losses  due  to  air  movement  through  the  inside  of  the  cylinder. 
The  temperature  of  the  cylinder  was  measured  at  the  same  point  along  the 
cylinder  as  the  amplitude. 

Vibration 

The  frequency  of  vibration  was  set  on  the  dial  of  the  audio 
frequency  oscillator,  but  a  strobotac  was  used  to  measure  the  frequency 
of  vibration  on  the  cylinder  because  the  scale  of  the  strobotac  was 
marked  in  smaller  graduations  than  the  dial  of  the  oscillator. 

The  amplitude  of  vibration  was  measured  from  the  photograph  taken 
of  the  cylinder  motion  at  the  time  the  temperature  was  recorded. 

Preliminary  to  data  runs  for  each  cylinder  installation,  a  photograph 
was  made  of  the  stationary  test  cylinder.  This  photograph  was  used 
to  determine  the  scale  factor  for  that  specific  cylinder  and  optics 
adjustment  by  comparing  the  width  of  the  cylinder  on  the  photograph  to 
the  actual  diameter  of  the  cylinder  (see  Figure  6).  This  factor  was 
used  for  subsequent  test  runs  to  compute  the  true  amplitudes  of 
vibration.  Amplitudes  were  measured  from  the  photographs  with  an 
optical  measuring  instrument  which  was  graduated  to  0.0001  in. 

Electrical  Heating  Power 

The  electrical  power  for  heating  the  test  cylinder  was  adjusted 
by  varying  the  resistors  while  monitoring  the  cylinder  temperature 
readings.  When  the  power  was  finally  adjusted  to  give  the  desired 
cylinder  temperature,  the  power  was  obtained  from  the  current  and 
voltage  readings  using  the  formula  P  =  El. 
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IV.  Experimental  Procedure 
Testing  with  Cylinder  Stationary 

With  the  cylinder  stationary  and  the  optical  equipment  adjusted, 
the  cylinder  was  heated  to  a  predetermined  surface  temperature,  and  a 
photograph  was  taken  of  the  stationary  cylinder  for  determining  the 
scale  factor.  The  temperatures  and  electrical  power  measurements  were 
recorded  for  determining  the  free  convection  heat  transfer  rate. 

Testing  with  Cylinder  Vibrating 

The  cylinder  was  vibrated  at  a  selected  frequency  and  amplitude 
while  the  surface  temperature  of  the  cylinder  was  maintained  at  the 
predetermined  value  by  adjusting  the  electrical  power  input.  Since  the 
surface  temperature  of  the  cylinder  responded  quickly  to  changes  in 
power,  the  temperature  usually  stabilized  quickly  as  was  observed  on 
the  recorder,.  When  the  temperature  was  stabilized,  the  measurements 
of  temperature ,  electrical  power,  frequency  and  amplitude  of  vibration, 
and  barometric  pressure  were  recorded. 

Visualization  of  Boundary  Layer 

Special  test  runs  using  the  0.25  in  diameter  cylinder  heated  to 
200  F  above  room  temperature  were  performed  at  vibration  intensities  in 
the  vicinity  of  the  critical  Reynolds  number  to  photograph  the  behavior 
of  the  boundary  layer.  The  same  procedure  as  in  stationary  and-  vibratory 
testing  wac  followed  with  one  exception — the  Schlieren  optical  system 
was  used.  When  the  desired  image  was  obtained,  the  test  was  conducted 
as  usual,  and  a  photograph  of  the  boundary  layer  was  taken  at  the  time 
the  variables  weie  recorded. 
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V.  Results 

The  results  of  this  investigation  are  presented  in  both  tabular 
and  graphical  form.  Measured  variables  and  calculated  parameters  for 
all  test  runs  are  contained  in  Tables  I  through  III  in  Appendix  E. 

The  calculated  parameters  are  presented  graphically  in  Figures  7 
through  9. 

There  exists  for  each  cylinder  diameter  a  critical  Reynolds 
number.  Vibration  intensities  below  this  Reynolds  number  give  heat 
transfer  rat^es  expressed  as  Nusselt  numbers  that  are  obtained  in  a 
free  convection  process.  Vibration  intensities  above  this  Reynolds 
number  give  Nusselt  numbers  that  are  obtained  in  a  forced  convection 
process.  The  Nusselt  numbers  corresponding  to  Reynolds  numbers  above 
the  critical  Reynolds  number  were  generally  along  the  forced  convection 
curve  of  McAdams  to  a  Reynolds  number  of  1100  which  was  the  maximum 
obtained  in  this  investigation.  The  critical  Reynolds  number  was  either 
on  or  near  the  conventional  forced  convection  curve. 

The  heat  transfer  rates  expressed  in  Nusselt  numbers  at  a  given 
Reynolds  number  were  not  different  for  a  cylinder  surface  temperature 
of  100  F  and  200  F  above  room  temperature. 

The  heat  transfer  rates  expressed  as  Nusselt  numbers  showed  no 
consistent  trend  for  the  effect  of  cylinder  diameter.  The  rates  were 
not  independent  of  cylinder  diameter,  nor  did  they  consistently  increase 
or  decrease  with  cylinder  diameter.  Rather,  the  results  indicate  that 
the-  rate  increases  with  cylinder  diameter  up  to  a  certain  diameter  and 
then  decrease  with  larger  cylinder  diameters. 

The  heat  transfer  rate  is  independent  of  the  test  position  along 


12 


6AW/ME/65-3 


the  cylinder,  and  is  only  dependent  upon  the  vibration  intensity  at 
the  test  point. 

The  boundary  layer  gradually  became  thinner  and  showed  an  increase 
in  turbulence  above  the  cylinder  as  the  vibration  intensity  was 


increased. 
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VI .  Discussion  of  Results 

The  variation  of  the  heat  transfer  rate  with  vibration  intensity 

•  Nvi 

is  presented  for  each  cylinder  in  the  non-dimensional  forms  - -  and  Re 

Pr  ♦ 3 

in  Figures  7,  8,  and  9.  Vibration  intensity  is  expressed  in  Reynolds 
numbers,  and  the  highest  Reynolds  number  shown  for  each  cylinder  was 
the  maximum  possible  with  the  test  apparatus  used. 

Since  the  primary  objective  of  this  investigation  was  to  test  to 
the  highest  possible  Reynolds  numbers,  most  of  the  data  were  obtained 
with  the  cylinders  vibrated  in  the  fundamental  mode,  for  this  mode 
produced  the  highest  Reynolds  numbers.  Some  data  were  obtained  with 
the  two  smaller  cylinders  vibrating  in  their  second  and  third  harmonics. 
The  0.25  in  diameter  cylinder  could  be  vibrated  only  in  the  fundamental 
mode.  The  heat  transfer  rate  from  each  individual  cylinder  was 
independent  of  all  modes,  positions  along  the  cylinder,  and  cylinder 
temperatures.  It  depended  only  upon  the  vibration  intensity  but  showed 
some  variation  between  the  different  diameter  cylinders.  The  broken 
lines  on  the  graphs  show  the  limits  of  scatter  in  the  data  from  which 
the  conventional  forced  convection  curve  of  McAdams  is  derived.  Most 
of  the  data  of  this  investigation  fall  within  the  limits  of  scatter 
shown  by  these  lines. 

The  critical  Reynolds  number  for  each  cylinder  falls  either  on 
the  conventional  forced  convection  curve  or  near  it.  Beyond  the 
critical  Reynolds  number,  the  heat  transfer  rate  increases  with  Reynolds 
number  and  generally  parallels  the  forced  convection  curve  to  a  Reynolds 
number  of  at  least  1100  which  was  the  maximum  obtained. 

An  analysis  of  the  data  for  the  three  cylinders  shows  no  consistent 
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trend  of  the  effect  of  cylinder  diameter  on  the  heat  transfer  rate 
expressed  as  a  Nusseit  number.  The  Nusselt  number  at  a  given  Reynolds 
number  for  the  0.120  in  diameter  cylinder  was  significantly  larger  than 
for  the  0.072  in  diameter  cylinder,  but  the  Nusselt  number  for  the  0.25 
in  diameter  cylinder  was  slightly  less  than,  for  the  0.120  in  cylinder. 

Comparison  with  Other  Authors 

A  graphical  comparison  with  the  results  of  other  authors  who 

investigated  heat  transfer  from  vibrating  cylinders  in  free  convection 

Nu 

is  presented  in  the  non-dimensional  forms  - -  and  Re  in  Figure  12.  An 

Pr*  3 

average  curve  is  shown  for  the  data  of  this  investigation  for  purposes 
of  comparison,  and  a  discussion  is  presented  in  the  following  paragraphs 
comparing  this  curve  to  the  data  curve  of  other  authors. 

McAdams .  The  data  for  the  0.072  in  diameter  cylinder  show  a 
heat  transfer  rate  which  is  an  average  of  8.6%  below  the  conventional 
forced  convection  curve  but  generally  parallel  to  the  curve.  Transition 
from  the  free  convection  heat  transfer  rate  to  the  forced  rate  occurs 
at  a  Reynolds  number  of  25;  whereas,  the  transition  on  the  conventional 
forced  convection  curve  occurs  at  a  Reynolds  number  of  7. 

The  data  for  the  0.120  in  diameter  cylinder  are  within  a  few  percent 
of  the  conventional  forced  convection  curve  and  generally  parallels  it. 

The  Reynolds  number  at  which  the  heat  transfer  rate  changes  from  a  free 
convection  rate  to  a  forced  convection  rate  appears  to  coincide  with  the 
Reynolds  number  at  which  the  free  convection  rate  crosses  the  conventional 
forced  convection  curve. 

The  data  for  the  0.25  in  diameter  cylinder  fall  along  the  con¬ 
ventional  forced  convection  curve  except  at  Reynolds  numbers  near  1000 
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where  the  heat  transfer  rates  are  approximately  10%  higher  than  the 
conventional  curve.  The  heat  transfer  rate  changes  from  that  obtained 
in  a  free  convection  process  to  that  obtained  in  a  forced  convection 
process  at  a  Reynolds  number  that  nearly  coincides  with  the  Reynolds 
number  at  which  the  free  convection  rate  crosses  the  conventional  forced 
convection  curve. 

The  data  as  a  whole  fall  within  a  band  which  is  approximately  10% 
either  side  of  the  conventional  forced  convection  curve  of  McAdams, 
and  which  is  generally  parallel  to  this  forced  convection  curve. 

Shine  and  Jarvis.  A  comparison  of  the  data  for  the  0.072  in 
diameter  cylinders  shows  the  heat  transfer  rates  for  the  cylinder  of 
this  investigation  to  be  26%  higher  than  the  rates  obtained  by  „Shine 
and  Jarvis,  though  the  slopes  of  the  two  curves  are  identical.  A 
similar  difference  occurred  between  the  results  of  Neely  and  the 
results  of  Shine  and  Jarvis.  As  did  Neely,  this  author  examined  the 
data,  equations,  and  computations  of  Shine  and  Jarvis  but  was  unable 
to  find  an  explanation  for  the  difference. 

Russ  and  Neely.  Since  Russ  and  Neely  both  used  similar  test  equip¬ 
ments  and  Neely's  investigation  was  in  extension  of  the  investigation 
performed  by  Russ,  a  comparison  is  made  jointly  with  the  results  of  these 
authors.  They  vibrated  the  entire  length  of  the  cylinder  uniformly  and 
transverse  to  the  length.  Russ  obtained  a  maximum  Reynolds  number  of 
152  using  cylinders  of  diameters  0.085  in,  0.25  in,  and  0.75  in  which 
were  heated  from  125  to  167  F.  Neely  obtained  a  maximum  Reynolds  number 
of  896  using  cylinders  of  diameters  0,085  in  and  0.25  in  which  were 
heated  from  60  to  120  F  above  room  temperature  (Ref  6:6). 

The  only  difference  between  the  data  of  Russ  and  Neely  and  the 
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data  of  this  investigation  is  in  the  range  of  Reynolds  numbers  from 
0  to  100.  The  data  of  this  investigation  show  that  the  heat  transfer 
rate  changes  from  that  obtained  in  a  free  convection  process  to  that 
obtained  in  a  forced  convection  process  at  the  conventional  forced 
convection  curve  and  generally  parallels  this  curve;  whereas,  the  heat 
transfer  rate  from  the  investigations  of  Russ  and  Neely  changes  at  a 
60%  higher  Reynolds  number,  and  the  transition  is  abrupt.  Also,  in  this 
range  of  Reynolds  numbers,  the  heat  transfer  rates  of  Russ  and  Neely 
are  as  much  as  50%  lower  than  the  rates  of  this  investigation.  Russ 
and  Neely  computed  the  free  convection  heat  transfer  rates  indirectly; 
whereas,  the  rates  were  directly  measured  in  this  investigation. 

Further,  Schlieren  photographs  taken  of  the  boundary,  layer  in  the 
vicinity  of  the  critical  Reynolds  number  show  ai  gradual  transition  in 
the  boundary  layer  which  correlates  with  the  data  of  this  investigation 
(see  Figure  13). 

Above  a  Reynolds  number  of  100  the  heat  transfer  rates  are 
identical  for  the  0.085  in  and  0.25  in  diameter  cylinders  of  Neely 
compared  to  the  0.120  in  and  0.25  in  diameter  cylinders  respectively 
of  this  investigation  (see  Figures  10  and  11). 

Fand  and  Kaye.  Fand  and  Kaye  used  a  0.875  in  diameter  cylinder  for 
which  they  obtained  Reynolds  numbers  from  955  to  2190  with  the  surface 
temperatures  ranging  from  25  to  185  F  above  room  temperature  (Ref  6:5). 
The  range  of  Reynolds  numbers  obtained  during  the.  investigation  by  Fand 
and  Kaye  began  approximately  at  the  maximum  value  pf  Reynolds  number 
obtained  during  this  investigation.  Though  the  heat  transfer  rates  at 
this  Reynolds  number  differ  between  the  two  curves  by  16%,  the  slopes  of 
the  curves  are  parallel  to  the  conventional  forced  convection  curve. 
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The  16%  difference  between  the  curves  possibly  indicates  an  effect  of 
cylinder  diameter  on  the  heat  transfer  rate.  The  data  from  this  inves¬ 
tigation  show  the  rate  increasing  from  the  0.072  in  diameter  cylinder 
to  the  0.120  in  diameter  cylinder,  and  then  decreasing  slightly  between 
the  0.120  in  diameter  cylinder  and  the  0,25  in  diameter  cylinder.  The 
results  for  the  0.875  in  diameter  cylinder  presented  by  Fand  and  Kaye 
show  a  further  decrease  in  the  heat  transfer  rate  between  the  0.25  in 
diameter  cylinder  and  the  0.875  in  diameter  cylinder.  Therefore,  the 
results  of  this  investigation  in  conjunction  with  the  results  of  Fand 
and  Kaye  point  to  a  possible  trend  of  the  effect  of  cylinder  diameter 
on  the  heat  transfer  rate.  This  trend  would  indicate  that  there  is  a 
certain  diameter  cylinder  between  0.072  in  and  0.25  in  diameter  that 
gives  a  maximum  heat  transfer  rate  for  all  Reynolds  numbers. 

Overall  Correlation  with  Other  Authors .  This  investigation 
1)  provides  what  may  be  a  more  accurate  behavior  of  the  heat  transfer 
rate  in  the  vicinity  of  the  critical  Reynolds  number,  2)  substantiates 
the  heat  transfer  rates  in  the  range  of  Reynolds  numbers  from  100  to 
1000,  and  3)  further  provides  substantiation  of  previous  results  above 
a  Reynolds  number  of  1000  because  of  close  correlation  with  these 
previous  results.  Therefore,  the  results  of  this  investigation  in  con¬ 
junction  with  the  results  of  Russ,  Neely,  Shine  and  Jarvis,  and  Fand 
and  Kaye  show  that  the  heat  transfer  rate  from  a  cylinder  of  diameter 
from  0.072  to  0.875  in  vibrated  in  air,  in  any  way  transverse  to  the 
length,  at  any  temperature  up  to  approximately  200  F  above  ambient 
temperature  will  change  from  being  a  free  convection  process  to  a 
forced  convection  process  and  follow  the  well-established  forced 
convection  curve  within  ±10%  up  to  a  Reynolds  number  of  at  least  2190. 


18 


GAW/ME/65-3 


Accuracy  of  Results 

Because  the  magnitude  of  the  heat  transfer  rates  and  vibration 
intensities  were  generally  much  larger  than  the  errors  associated  with 
the  variables,  the  errors  were  relatively  insignificant.  Certainly * 
the  errors  incurred  in  the  experiment  were  no  more  significant  in 
calculating  the  dimensionless  parameters  than  the  errors  in  reading  the 
slide  rule. 

Errors  which  might  have  been  caused  by  room  convective  air 
currents  and  conductive  heat  losses  from  the  cylinder  were  avoided  by 
enclosing  the  test  assembly  in  a  suitable  hood  and’  by  insulating  the 
cylinder  from  all  objects  except  the  copper  electrical  power  connections. 
The  cylinder  was  uniformly  heated  and  the  longitudinal  temperature 
gradient  was  nearly  zero  along  the  iength  of  the  cylinder  within  4  in 
either  side  of  the  test  point. 

The  maximum  error  which  occurred  in  the  Nusselt  number  based  on  all 
combined  errors  was  3.7%  for  the  0.072  in  cylinder  at  the  lowest 
vibration  intensity.  Likewise  the  maximum  error  in  the  corresponding 
Reynolds  number  was  2.0%.  The  errors  decrease  with  increasing  vibration 
intensities,  cylinder  diameter,  and  temperature  difference.  A  detailed 
discussion  of  errors  is  presented  in  Appendix  D. 

Test  Limitations 

The  scope  of  this  investigation  was  limited  by  cylinder  vibration 
characteristics,  cylinder  failure  at  higher  temperatures,  and  the 
capability  of  the  mounting  assembly  to  support  a  more  massive  cylinder. 

Higher  vibration  intensities  within  a  vertical  plane  were  not 
obtained  because  the  cylinder  motion  changed  from  two-dimensional  to 
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three-dimensional  beyond  a  critical  intensity.  Since  this  investigation 
was  restricted  to  motion  in  a  vertical  plane,  the  vibration  intensities 
achieved  were  limited  by  the  intensities  at  which  motion  of  the  cylinder 
left  the  vertical  plane. 

Tests  to  attain  high  Reynolds  numbers  with  cylinder  surface  tempera¬ 
tures  above  290  F  were  not  possible  because  of  failures  of  the  0.25  in 
cylinder.  The  failures  occurred  too  frequently  at  higher  temperatures 
to  continue  testing. 

The  mass  of  the  mounting  assembly  was  not  adequate  to  support  a 
larger  diameter  test  cylinder  than  the  0.25  in  cylinder.  High  intensity 
vibration  of  the  0.25  in  cylinder  excited  the  mounting  assembly  which 
in  turn  caused  the  motion  of  the  cylinder  to  be  too>  unsteady  to  allow 
accurate  data  to  be  recorded. 
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VII.  Conclusions 

Conclusions  drawn  from  the  results  of  this  investigation  are: 

1.  The  heat  transfer  rate  (expressed  in  a  Nusselt  number)  from  a 
vibrating  cylinder  in  air  changes  from  that  obtained  in  a  free 
convection  process  to  that  obtained  in  a  forced  convection 
process  at  a  critical  Reynolds  number.  This  Reynolds  number 
increases  with  cylinder  diameter  and  occurs  at  or  near  the 
conventional  forced  convection  curve  of  McAdams,  and  beyond 
this  Reynolds  number,  the  heat  transfer  rate  will  follow  the 
conventional  forced  convection  curve  within  approximately 
~10%  up  to  a  Reynolds  number  of  at  least  1100.  Further, 
based  on  close  correlation  with  the  results  or  Fand  and  Kaye, 
the  heat  transfer  rate  can  be  expected  to  follow  the  conven¬ 
tional  forced  convection  curve  up  to  a  Reynolds  number  of  at 
least  2190. 

2.  The  heat  transfer  rates  expressed  in  Nusselt  numbers  at  a 
given  Reynolds  number  were  not  different  for  a  cylinder 
temperature  of  100  F  and  200  F  above  room  temperature. 

3.  The  heat  transfer  rate  is  independent  of  the  position  of  the 
test  point  along  the  cylinder.  It  is  only  dependent  upon 
the  vibration  intensity  at  the  point. 

4.  No  general  effect  of  cylinder  diameter  on  the  heat  transfer 
rate  was  determined  by  this  investigation;  however,  in 
conjunction  with  the  results  of  Fand  and  Kaye,  there  appears 

to  be  a  trend  for  the  rate  to  increase  up  to  a  certain  cylinder 
diameter  and  then  decrease  with  larger  diameter  cylinders. 
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Figure  1 

Schematic  Layout  of  Experimental  Equipment 
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1-Lgure  4 

Photograph  of  Hood  Which  Encloses  the 
Cylinder  and  Surrounding  Space 


Figure  5 

Photograph  of  Rotating  Mirror  on  Small  D.C. 
Electric  Motor  in  Front  of  Camera  Aperture 
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Figure  10 

Comparison  of  the  Heat  Transfer  Rates  with  Vibration  Intensity 
Between  the  0.120  in  Diameter  Cylinder  and  the  0.085  in 
Diameter  Cylinder  of  Neely's  investigation 
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Figure  11 

Comparison  of  the  Heat  Transfer  Rates  with  Vibration  Intensity 
Between  the  0.25  in  Diameter  Cylinder  and  the  0.25  in 
Diameter  Cylinder  of  Neely's  Investigation 
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Photographs  of  the  Boundary  Layer  on  Cylinder  While  Vibrating 
in  the  Vicinity  of  the  Critical  Reynolds  Number 


oft H/ntjf  00-0 


Appendix  B 


Equations 


Nusselt  Number 


Nu 


3.413 


IE  d 
A  k  AT 


B  watt 
watt  hr  ft2 


ft  hr  ft  F 
B  F 


Reynolds  Number 


ft  ft  hr  ft  lbm  3600  sec 
sec  lbm  12  ^.3  hr 


where  a  is  measured  in  inches. 


Re  = 


vdp  _  2afdp 
4  ~  V 


=  600 


afdp 
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Appendix  C 


Sample  Calculation 


The  following  calculations  are  based  on  run  #20  for  the  0.25  in 
diameter  cylinder.  The  corrected  values  of  the  measured  variables  are: 


Ta  =  79  F 

I  =  20.0  amps 

a  =  0.3740  in 

T w  =  179  F 

E  =  1.75  volts 

f  =  42.70  cps 

B.P.  =  29.385  in  Hg 


From  tables  in  Eckert  and  Drake  (Ref  1),  the  fluid  properties  were 
determined  as 

kf  =  0.16395  BTU/hr  ft  F 
\i  =  0.04777  Ibm/hr  ft 
Pr  =  0.7042  dimensionless 


The  following  computations  were  performed: 

Tw  +  Ta  179  +  7g 
Tf  =  =  129  F 


AT  =  Tw  -  Ta  =  179  -  79  =  100  F 


Qt  =  IE  =  (20.0)(1.78)  =  35.0  watts 


The  radiation  heat  loss  computed  from  the  following  formula 
assuming  the  hood  was  very  large  compared  to  the  diameter  of  the 
cylinder  (Ref  5:390). 
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,  TK.4  .  T,  h 

Qr  =  0.173  £  Aw  (jjj)  -  yi) 


■  °-173  *  °-65  *<**"§*  »•»>  Q"  -  (ii)" 


=  1.123  x  ( . 2055) (10. 25  -  0.039) 


=  .01465 


=  0.05  watt 


The  conduction  losses  were  assumed  to  be  zero  because  the  cross-section 
area  of  the  cylinder  was  only  0.0375  in2  and  the  temperature  gradient 
along  the  cylinder  was  too  small  to  be  significant.  Bednorz  showed 
experimentally  that  conduction  losses  are  insignificant  (Ref  1:12). 

The  convection  heat  transfer  was 


Qc  =  Qt  -  Qr  =  (35.00  -  0.05)  =  34.95 


Nusselt  Number 


Nu  = 


0.333  Qt 
k  AT 


0.333  x  35.Q 
0.016395  x  100 


=  7.11 


Nu  _  7.11 
pr.3  "  .902 


=  7.88 


Reynolds  Number 


12.5  af  - 
VI 


=  12.5  X  0.374  >  42.7  x 


0.0661 

0.04777 
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The  overall  error  in  the  Nusselt  n camber  (see  Appendix  D) 


ANu  = 


3.413  d 
A  k 


A  (AT.) 
(AT)2 


The  overall  err>or  in  the  Reynolds  number 
Re  =  600  —  (fpAa  +  apAf  +  afAp) 

Aa  =  0.002  Af  =  0.617  Ap  *  0  p  =  .0661 

0  25 

ARe  =  600  -i^-x-oV 64777  (0‘002  x  42-7  x  -0661  +  0.374  x  0-617  x  0.0661) 
=  261  x  13.0209 


=  5.45 


The  maximum  overall  error  will  be  different  for  each  run  and  will 
decrease  as  the  size  of  the  cylinder  increases,  and  as  the  magnitude 
of  the  Nusselt  number  and  Reynolds  number  increases  the  percent  error 
decreases. 
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Appendix  D 


Analysis  of  Errors 

Both  human  and  instrument  errors  were  possibly  introduced  into  the 
raw  data  during  the  measurement  of  each  of  the  seven  variables.  An 
analysis  of  the  individual  errors  is  presented  in  the  following 
paragraphs . 

Power  Measurement 

Electrical  power  input  to  the  cylinder  was  computed  using  the 
formula  P  a  El. 

Though,  the  maximum  error  of  the  ammeter  and  voltmeter  was  1%  of 
the  full  scale  reading,  the  effect  of  this  error  was  reduced  by 
calibrating  the  instruments  against  a  meter  which  had  a  maximum  error 
of  1/4  of  1%  of  the  full  scale  reading.  The  only  other  error  was  the 
human  error  in  reading  the  scales  of  the  instruments.  The  accuracy  of 
data  read  on  the  instruments  are  presented  in  Table  D-l. 

Vibration  Intensity 

Vibration  intensity  is  the  product  of  amplitude  of  vibration  and 
frequency  of  vibration,  and  errors  were  possible  in  measuring  each  of 
these  factors. 

Frequency.  The  frequency  of  vibration  was  measured  with  a  Strobotac 
which  was  calibrated  prior  to  each  series  of  test  runs.  When  calibrated, 
the  accuracy  of  the  instrument  was  - 1 %  of  the  high  scale  which  amounts 
to  a  possible  error  of  ±0.617  cylces  per  second.  The  total  error  in 
frequency  data  is  shown  in  Table  D-i. 
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Amplitude ♦  The  amplitude  of  vibration  was  measured  on  photo¬ 
graphic  film.  The  possible  contributions  to  this  error  are  errors  in 
measuring  the  amplitude  of  the  cylinder  motion  from  the  photographs. 

The  amplitudes  were  scribed  with  a  knife  and  read  on.  a  micro  compensator 
which  can  be  read  to  the  fourth  decimal  place.  Thus  the  primary  error 
in  measuring  from  the  photograph  was  in  scribing  the  amplitude  on  the 
photograph  which  contributed  an  error  of  ±0.001  in.  The  optical 
system  was  calibrated  to  within  an  accuracy  of  ±0.001  in  error. 

Therefore,  the  maximum  total  error  in  obtaining  the  amplitude  was 
±0.002  in.  Hence,,  the  total  error  introduced  in  obtaining  the  vibra¬ 
tion  intensity  was  ±0.89  in/sec  as  computed  from  A(af)=  =  Aaf  +  aAf. 

Temperature 

The  temperature  measuring  system,  iron-constantan  thermocouples 
and  a  Bristol's  Dynamaster  Recorder  potentiometer,  was  calibrated 
against  a  mercury- in-glass  thermometer  in  an  ice  bath,  in  ambient  air 
(75  F),  and  in  boiling  water.  The  65  to  212  F  scale  of  the  recorder 
used  during  the  tests  with  a  AT  of  100  F  agreed  with  the  thermometer 
as  close  as  the  two  scales  could  be  read.  The  40  to  320  F  scale  used 
during  the  test  with  a  AT  of  200  F  agreed  with  the  thermometer  at  32  F 
and  75  F  but  differed  at  212  F  by  2  F.  This  error  was  applied  to 
subsequent  readings  above  212  F.  The  temperature  response  to  power 
changes  and  vibration  intensity  was  immediate.  This  can  be  attributed 
to  a  positive  contact  which  was  maintained  between  the  cylinder  surface, 
a  small  size  thermocouple  junction,  and  the  small  heat  capacity  of  the 
cylinder  and  thermocouple.  Because  of  the  quick  response  in  temperature, 
the  temperatures  read  were  well  stabilized,  thus  minimizing  any  error 
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in  readings.  The  stability  of  cylinder  and  room  temperatures  was 
further  increased  by  the  hood  which  enclosed  the  space  surrounding  the 
test  cylinder.  Therefore,  a  close  estimate  of  the  accuracy  with  which 
temperatures  were  obtained  is  1/2  F — the  closest  the  scale  could  be 
read. 


Overall  Error 

The  error  in  Nusselt  number  due  to  errors  in  variables  is 


ANu 


3.41  d 
A  k  ' 


AE 

AT 


+  IE 


A  (AT)' 
‘(AT.)"2 


The  error  in  Reynolds  number  due  to  errors  in  variables  is 

r\ 

ARe  =  600  —  (Aafp  +  aAfp  +  afAp) 

The  maximum  overall  error  computed  for  the  test  run  which  should 
have  the  largest  error  during  this  investigation  was  a  3.72%  error  in 
Nu  and  a  2.0%  error  in  Re. 


Appendix  E 


Experimental  Data 
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Personally  Identifiable 
Information  Redacted 


He  was  graduated  from  the  University  of  Alabama  in  1952  and  was 
commissioned  as  a  second  lieutenant  in  the  United  States  Air  Force 
through  the  ROTC  program.  At  that  time  he  entered  pilot  training,  and 
upon  completion  of  training  was  assigned  duty  with  the  44th  Fighter 
Bomber  Squadron  at  Clark  AFB,  P.  I.  Upon  completing  an  eighteen  month 
tour  of  duty,  with  the  44th  FBS,  he  was  assigned  to  duty  with  the  Air 
Force  Systems  Command  at  Elgin  AFB,  Florida  and  Wright -Patterson  AFB, 
Ohio.  Subsequently,,  he  entered  resident  training  for  Graduate 
Aeronautical  Engineering  (Air  Weapons  Option)  at  Wright-Patterson  AFB, 
Ohio. 
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